Astronomy & Astrophysics manuscript no. nymarkXraySN 


February 5, 2008 


(DOI: will be inserted by hand later) 





X-ray emission from radiative shocks in Type II supernovae 

Tanja K. Nymark, Claes Fransson, and Cecilia Kozma 

Stockholm Observatory, AlbaNova University Center, SE-106 91 Stockholm, Sweden 
Received 8 Sep 2005/ Accepted 21 Oct 2005 

Abstract. The X-ray emission from the circumstellar interaction in Type II supernovae with a dense circumstellar 
medium is calculated. In Type IIL and Type Iln supernovae mass loss rates are generally high enough for the 
region behind the reverse shock to be radiative, producing strong radiation, particularly in X-rays. We present 
a model for the emission from the cooling region in the case of a radiative reverse shock. Under the assumption 
of a stationary flow, a hydrodynamic model is combined with time dependent ionization balance and multilevel 
calculations. The applicability of the steady state approximation is discussed for various values of the ejecta density 
gradient and different sets of chemical composition. We show how the emerging spectrum depends strongly on 
the reverse shock velocity and the composition of the shocked gas. We discuss differences between a spectrum 
produced by this model and a single-temperature spectrum. Large differences for especially the line emission are 
found, which seriously can affect abundance estimates. We also illustrate the effects of absorption in the cool 
shocked ejecta. The applicability of our model for various types of supernovae is discussed. 
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' 1. Introduction 

. The emission at late times is for most Type IIL (Type 
' II linear) and Type Iln (Type II narrow line) super- 
. novae dominated by the interaction between the super- 
' nova ejecta and the circumstellar medium. The collision 
between the two creates two shocks; one moving out- 
ward into the circumstellar medium and a reverse shock, 
which moves backward into the ejecta (Fig. Because 
' of the high temperatures involved, most of the emis- 
. sion from the interaction region is radiated as X-rays. 
' The strong emission heats and ionizes the circumstellar 
medium and ejecta, leading to line emission in the optical 
and UV fFransson ll984l Chevalier & Fransson 1994, here- 
after [^^^. This can then act as probes of the ejecta and 
the shock. Another characteristic feature of the interac- 
tion is strong radio emission caused by synchrotron emis- 
sion behind the outer shock (e.g., Sramek et al. I2003|l . A 
review of circumstellar interaction in supernovae is given 
by Chevalier & Fransson (2003 ). 

Chevalier H1982al[l982b|l proposed circumstellar inter- 
action as an explanation of the observed radio light curves 
of SN 1979C and SN 1980K, and used a self-similar for- 
mulation to analyze the effects of interaction between su- 
pernova ejecta and the ambient medium. Different mech- 
anisms for production of X-rays were also discussed. In 
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Fransson ri984P the X-ray emission from the reverse and 
circumstellar shocks was discussed in some detail. In par- 
ticular, it was argued that in most cases the reverse shock 
would be radiative, a fact we will exploit in this paper. 
This also has the consequence that a cool, dense shell 
will form between the two shocks, which is important for 
both the X-ray and optical emission. A detailed study of 
the effects of circumstellar interaction in Type II super- 
novae was done by Chevalier & Fransson l|CF94|l . while 
the specific case of SN 1993J was investigated by Fransson 
et al. (1996, hereafter iFLC96|l . as well as by Suzuki & 
Nomoto (1995). 

Radiative shocks in supernova remnants were first 
modeled by Cox l|1972|) . Later, improved models were 
made by Raymond H1979|l . Dopita I|1976[I1977|) and others. 
Because of the lower densities, the shocks are, however, 
only radiative for velocities ^ 200 km s~^, substantially 
lower than the velocities studied in the present paper. High 
density models, suitable for active galactic nuclei (AGN) 
and young supernova remnants were made by Plewa &; 
Rozyczka (fUMIl and Plewa lfT!Ml [TM^ . These models 
are, however, mainly concerned with the hydrodynamics, 
with a simplified treatment of the emergent spectra. 

The goal of the work in this paper has therefore been 
to develop a model for the X-ray emission from the reverse 
shock. This combines the hydrodynamic structure with a 
spectral code which, by including updated atomic data 
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and multi-level solutions to most ions, can model the line 
emission in greater detail than has been done before. 

In Sect. 121 we give a brief summary of the most im- 
portant observations of circumstellar interaction, and in 
particular X-ray observations. Our model is described in 
Sects. |31through|Hl while Sect.[7|describes the results of our 
calculations, including the contribution from the circum- 
stellar shock. The conclusions are summarized in Sect.|H| 
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Fig. 1. Schematic picture of the different regions in the 
interaction of supernova ejecta with the progenitor wind. 
The figure is not to scale. The contact discontinuity be- 
tween the shocked wind and shocked ejecta is Rayleigh- 
Taylor unstable. 



2. X-ray emission from supernovae 

The main observational indications of circumstellar inter- 
action are the detection of narrow optical and UV emis- 
sion lines from ionized circumstellar gas, broad, box-like 
emission lines from ejecta ionized by radiation from the 
interaction region, and a wavelength dependent turn-on 
of the observed radio emission. The latter is explained as 
a result of either free-free absorption or synchrotron self- 
absorption. Another effect of circumstellar interaction is 
that the supernova at late epochs displays a flat optical 
light curve, long after the emission caused by radioactive 
decay has faded away. For Type Iln supernovae circum- 
stellar interaction dominates the optical spectrum already 
from early epochs. We refer to Filippenko H1997|l for a re- 
view of the optical observations. 

Besides the radio emission, the X-rays show the clear- 
est evidence for circumstellar interaction. Because the X- 



ray emission is the focus of this paper, we will briefly re- 
view the observational evidence for this. 

Until now, ~ 24 supernovae have been detected in X- 
rays, and for 10 of these the data are of sufficiently 
good quality to be studied in detail. We will concentrate 
on those for which good data exist, and circumstellar in- 
teraction is thought to be important. Comprehensive re- 
views of the existing X-ray observations of supernovae can 
be found in Schlegel H1995|l and Immler & Lewin [20^^. 

The first supernova to be detected in X-rays was 
SN 1980K, which was observed by the Einstein satellite 35 
days after outburst (Canizares et al. ll982|l . At the time of 
the next observation at 82 days, the emission had decayed 
by at least a factor of two. The emission was found to be 
rather soft, with kT ^ 0.5 keV, but no spectra were ob- 
tained. Also SN 1979C was observed with Einstein at early 
epochs, but no source was seen (Panagia et al. .I98IHl . After 

16 years it was, however, detected as a strong source 
by ROSAT (Immler et al. HMSjl . It was observed by the 
ASCA satehite in 1997 (Ray et aLEHDl) and by Chandra 
in 1999 and 2000 (Kaaret'^UDT), with a roughly constant 
X-ray fiux in the period of 16-20 years after outburst. The 
emission was found to be soft, and likely to originate in 
the shocked ejecta. The reason for the absence of strong 
X-rays at early epochs is most likely absorption by the 
cool, dense shell formed by the reverse shock l|CF94|l . 

The Type Iln SN 1986J in NGC 891 was observed with 
ROSAT and ASCA from nine to thirteen years after out- 
burst (Houck et al. 1998). During this period it was found 
to have a fast decline of the X-ray emission, with oc t~^. 
From the ASCA observations it was found to have a rather 
hard spectrum {kT ~ 5 — 7 keV) with a clear Fe K emission 
line at 6.7 keV, with a width < 20 000 km s'^ 

The first supernova to be identified through its X-ray 
emission was the Type Iln SN 1978K. Strong radio and 
optical emission had been detected earlier (Ryder 119931 
Dopita & Ryder'1990'), but it was with the detection of X- 
rays by ROSAT in 1992 fRvder .1993) that it was identified 
as a supernova. Extensive observations in radio, optical, 
UV and X-rays were carried out in the late 1990s, and 
observations by Newton-XMM in October 2000, and by 
Chandra two years later, indicated a constant X-ray flux 
over the past decade (Schlegel et al. 12004(1 . In the early 
X-ray spectra no emission lines were found, but in the 
Chandra data lines from Si are evident. 

The best observed case of circumstellar interaction is 
SN 1993J in M81. Since its discovery in March 1993, it 
has been extensively studied in all wavebands. SN 1993 J 
showed signs of circumstellar interaction early on, first 
through the detection of radio emission on day 5 (Weiler 
et al. 11993(1 . and shortly thereafter through X-ray obser- 
vations. The first X-ray detection was made by ROSAT as 
early as 6 days after discovery (Zimmerman et al. I1994(l . 
and the first observation by ASCA came two days later. 
The early ASCA observations showed a strong Fe K line at 

^ A complete list of the X-ray supernovae may be found at 
http://lheawww.gsfc.nasa.gov/users/immler/supernovae_list. 
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6.7 keV, with a width of 0.3 ±0.3 keV, which faded within 
weeks after the explosion (Kohmura et al. 119941 Uno et 
al.^UUTi. Both ROSAT and ASCA observations showed a 
hard spectrum during the first two months. This was con- 
firmed by OSSE observations on Compton/GRO, where 
hard X-ray emission with kT ~ 50 — 150 keV, was de- 
tected in two observations, 9.9-15.4 days after explosion 
and again between days 23.7-36.9 (Leising et al. |1994). 
Three months after explosion it was too weak to be ob- 
served by OSSE. When it was re-observed by ASCA and 
ROSAT after six months a soft component was found to 
dominate the spectrum. This transition from a hard to a 
soft spectrum was explained by FLC96| as a result of de- 
creasing optical depth in the cool shell between the two 
shocks, which allows the soft X-rays from the reverse shock 
to escape. 

Very late Chandra observations on day 2594 yielded 
a detailed spectrum of SN 1993J in the 0.3-8.0 keV band 
(Swartz et al. .2003j . This spectrum showed evidence of 
Fe L emission, as well as emission from H-like Mg and Si 
and He-like Mg, Si and S. This is taken by Swartz et al. to 
indicate the presence of two low-temperature components, 
while the high energy (<;1.5 keV) continuum emission is 
thought to point to a high-temperature component as well. 

The Type Iln SN 1995N was detected by ROSAT in 
1996. The explosion date is unknown, but Fransson et 
al. H2002|l argue that based on the spectrum a likely ex- 
plosion date was July 1994. It was observed by ROSAT 
in 1996 and 1997, by ASCA in 1998 (Fox et al. IMmI 
and by Chandra in 2004 (Chandra et al. "2DDF). The X- 
ray luminosity was found to decline between the first two 
observations. Fox et al. (2000) found that the luminosity 
increased again by a factor of two from 1997 to 1998, al- 
though a reanalysis by Chandra et al. (2005) of the ASCA 
observations indicate that the source was contaminated 
and that both the ASCA and Chandra observations are 
consistent with a linear decline of the X-ray luminosity. 
Chandra et al. 1)2005(1 also find clear evidence of Ne line 
emission around 1 keV. 

The other well observed Type Iln SN 1998S in 
NGC 3877 was discovered at an early phase, and has 
been well studied at all wavelengths. It was observed 
by Chandra on five occasions from day 678 to day 1324 
(Pooley et al. l2001|l . A number of emission lines from met- 
als were observed at energies above 1 keV, among them 
lines from Ne, Si, S and Fe. The strengths of the observed 
lines were claimed to indicate a high overabundance of 
the metals with respect to solar values. This would imply 
that the observed X-ray emission mainly comes from the 
shocked oxygen core. In this case the heavy elements must 
be mixed out to a high velocity, possibly caused by an as- 
pherical explosion (Pooley et al. 2001 j . The abundance 
analysis is, however, based on a one-temperature model of 
the spectrum. In Sect. ITTI we show that this assumption 
is highly questionable. 

SN 1999em was the first Type IIP supernova to be 
detected both in X-rays and radio. It was detected by 
Chandra three days after discovery, and was observed 



seven times in a period of almost two years. During this pe- 
riod the spectrum softened from a 5 keV spectrum to a 
~ 1 keV spectrum. Pooley et al. (|2001ll modeled the X-ray 
emission as resulting from interaction between supernova 
ejecta and the wind of a red supergiant, and found that 
the X-ray emission agrees well with this, and that both 
shocks were adiabatic. 

The most recent X-ray supernova is the Type IIP 
SN 2004dj, which was detected by Chandra in August 
2004, most likely two months after explosion. It was 
found to have a temperature of kT ^7 keV, and nearly 
three times the luminosity of SN 1999em (Pooley & 
Lewin I2004|l . Modeling of Type IIP supernovae is dis- 
cussed in detail by Chevalier et al. H2005|l . 

3. Hydrodynamics 

3.1. Progenitor wind 

The immediate surroundings of the exploding star consist 
of the stellar wind of the progenitor. The density of a 
wind with constant mass loss rate and velocity drops as 
r~^, but the wind parameters depend on the evolutionary 
state of the star. On the main sequence the star is a blue 
supergiant with a fast wind. It then evolves to become 
a red supergiant with a correspondingly slow and dense 
wind, with mass loss rate of M ~ 10~^ — 10~* Mq yi:~^ 
and a velocity of ^ 5 — 50 km s^^ (e.g., Salasnich 
et al. 11999(1 . Most massive stars explode in this stage as 
supernovae. Stars more massive than ~ 22 Mq (Meynet & 
Maeder 12003(1 may, however, evolve back to the blue and 
become Wolf-Rayet stars or, in some cases, end their lives 
as blue supergiants. 

As the star evolves, the wind will create an onion-like 
structure of circumstellar gas, with dense shells alternat- 
ing with low density shells. In some stages the progenitor 
may undergo pulsations (Heger et al. 11997(1 . which could 
result in a wind with a density differing from p oc r~^. The 
circumstellar medium can also be modified by mass loss in 
a binary system, as might have been the case for SN 1993J 
(e.g., Suzuki 19M: Suzuki & Nomoto ISMI Nomoto et 
al. 11995(1 . Because of this, the assumption that the den- 
sity drops as distance squared may not always hold. A 
more general description of the wind density is therefore 

= iJk iy) ' 

where i?, is a reference radius. For s = 2 this reduces to 
the expression for a steady wind. There are few cases for 
deviations from s — 2 winds (see Fransson & Bjornsson, 
120051 for a discussion). In some cases, the supernova may, 
however, interact with either clumpy winds or shell-like 
structures, for which s may differ from 2 (see Sect. I3.3)l . 

3.2. Ejecta 

The unshocked supernova ejecta are assumed to expand 
freely, so that the position of a specific mass element can 
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be expressed as R{t) = Vt, where the velocity V is con- 
stant. The structure of the ejecta reflects the structure of 
the progenitor and the explosion itself. 

A useful approximation to the density structure of the 
outer ejecta is 



Poj 



r 



(2) 



where rj is the power law index characterizing the outer 
ejecta. Analytical ejecta models by Matzner & McKee 
H1999|l indicate that rj = 11.9 for a red supergiant pro- 
genitor and r] = 10.1 for a blue supergiant progenitor. 
For this to be valid, the progenitor structure has to be 
close to a polytrope of index 3/2 (convective envelope) or 
3 (radiative envelope), respectively, which does not always 
hold. In particular, the models of Matzner & McKee do 
not apply to superadiabatic regions, which can be present 
in certain progenitors. This may have been the case for 
SN 1993J (Nomoto et al. HM^ . 

In more detailed hydrodynamical explosion models the 
power law index i] is usually in the range 7-12, but can 
be as high as ~ 20, as is the case in the ejecta models of 
SN 1993J. Evolutionary progenitor models (Shigeyama & 
Nomoto nMH Woosley & Weaver HMf^ indicate that a 
power law structure is a good approximation for the outer 
layers, and to some extent also in the inner parts, albeit 
with a different value for the power law index, 77. Roughly 
speaking, most one dimensional models show a steep inner 
power law, a shallow middle part and a steep outer part, 
although there is considerable structure, particularly in 
the inner parts. However, hydrodynamic instabilities at 
especially the H/He interface as well as the He/0 core 
can modify the structure appreciably (Miiller et al. 119911 
Iwamoto et aL HWI Kifonidis et al. lM)^ . 

During the first weeks or months, the reverse shock 
passes through the outermost steep part and enters the 
shallower gradient toward the middle. The shallow part 
is traversed quickly, and eventually the shock enters the 
inner steep part. As the shock moves into a region with 
a different density slope, the X-ray spectrum will change. 
The time at which this spectral change occurs, could help 
constrain the explosion models. We will discuss this in 
more detail in a later paper. 

3.3. Interaction region 

Chevalier fl982a) found a self-similar solution for the 
interaction between supernova ejecta and a stationary 
medium, which provided a complete solution to the physi- 
cal quantities of the two shocks (see also Nadvozhin ll985|l . 
The main parameters are the power law exponents, tj and 
s, of the ejecta and wind densities, the mass loss parame- 
ter M/vw and the shock velocity (or equivalently the total 
energy). 

These relations can be simplified further if one assumes 
that the interaction shell is very thin compared to its ra- 
dius (the thin shell approximation), which is usually an 



adequate approximation. This becomes exact in the limit 
of strong radiative cooling, which is the case for the reverse 
shocks discussed in this paper. The circumstellar shock is, 
however, in general adiabatic. A similarity solution de- 
scribing this case is discussed in ICF94I In the following 
the subscript '0' refers to quantities immediately behind 
the reverse shock, 'rev' refers to parameters of the shock 
itself, while the subscript 'cs' refers to the circumstellar 
shock. 

The density and temperature behind the reverse shock 
can in the thin shell approximation be written as 1|FLC96|I 



Tn 



(ry-3)(r;-4) 
(3-s)(4-s) 



Pc 



M (3 



Mcs iv - 3)' 



rTcs - 2.27 X lOV 



(3 



2 



K 



(3) 

(4) 



where /i is the mean mass per particle (i.e., electrons plus 
ions) behind the reverse shock, and fics is the correspond- 
ing value behind the circumstellar shock. The density be- 
hind the outer shock is related to the wind density by 
Pes — 4/ciw. Setting r — Rs = Vcjt in Eq. we can ex- 
press the post-shock ion number density as 



no = 1-2 X 10-^ 



; (77-3)(r;~4) A, 
(3-s)(4-s) paVM- 



(5) 



where /iA is the mean atomic weight, and we have defined 



A, =M^5/v^, = (M/10-5Mo)/(vw/10 km s'^). Note 
that A* — lOOA*, where A, was defined by Chevalier & 
Li (1999) for parameters typical of a Wolf-Rayet wind. Rg 
is the radius of the interaction region, which we here take 
to be the same as the position of the reverse shock, and Vcj 
is the maximum velocity of the ejecta, i.e., the velocity of 
the ejecta at the reverse shock front. Note that the radius 
of the circumstellar shock, i?cs, is a factor 1.2-1.3 larger 
than Rs 



iDF94i. 

The temperature behind the circumstellar shock is 



Tcs = 2.27 X lO^^c 



(^-3)^ 

'iv 



vi K, 



^2 "^4 



(6) 



where = V^j/lO^ km s~^. Values of p. and are given 
in Table n] for the different compositions considered in this 
paper. The composition of the circumstellar shock is the 
same as in the circumstellar medium, which we assume to 
be solar, with ~ 0.61. The mass swept up by the outer 
shock front is 



MR, 



(2-s) 



(3 - s)v^ 

and the mass swept up by the reverse shock is 



■n 



-Mr. 



The velocity of the reverse shock is 

_ 3-s 



(7) 



(8) 



(9) 
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which in combination with the shock jump condition gives 
the velocity of the post-shock gas 



(3~^) 
4(77 - s) 



(10) 



As an alternative to using a self-similar model, depen- 
dent on the assumption that the ejecta density is a power 
law, p ^ r~'', we have the option of using a more detailed, 
non-power law model for the ejecta as input. In this case 
we use the thin shell approximation to solve for the po- 
sition of the shock, and consequently the density, tem- 
perature and velocity of the ejecta at this position (see 
e.g. Chevalier & Fransson l2(JU3fl . The post-shock values 
are then found from the shock jump conditions. Once the 
state of the post-shock gas has been found, the treatment 
of the interaction is the same in both cases, as described 
below. 

An alternative scenario for the X-ray emission is shocks 
driven into dense circumstellar clumps (e.g. Chugai ri993|l . 
In this case the X-ray emission would originate in radiative 
shocks within each clump. The structure of the radiative 
shocks would, however, for the same shock velocity be the 
same as for the reverse shock considered here. An espe- 
cially interesting case of this type is the ring collision of 
SN 1987A (Michael et aJ. IMT^ . 

3.4. Reverse shock structure 

The structure of the shocked ejecta behind the reverse 
shock is found by solving the conservation equations. If 
we reduce the problem to one dimension, which is justified 
for a thin shock, and assume steady state, these equations 
reduce to 



— imv) = 
ax 



— I MATOu^iW + Pj 

ax 



= 



_d_ 

dx 



-^AninriiV H - 

2 7 — 1 



+ p ] V 



(11) 
(12) 

= A{T)nin„ (13) 



where rij is the ion number density and A is the cooling 
rate. These equations may be solved to give the structure 
of the radiative shell, once the values immediately behind 
the shock are known. For the thin shell model the initial 
density and temperature are given by Eqs. (PJ and Q), 
and the initial velocity by Eq. (|l()|l . The initial pressure is 
Po = (n-0,0 + n-o)A:To 

The model is specified by the density, temperature and 
composition behind the shock. The initial value of the 
density depends on which structure we assume for the 
ejecta (see Sect. 13. 3f) . Once the density structure is set, 
the pressure and velocity of each zone can be found from 
the continuity equations 



V = novon^ 



(14) 



p = /iAnii,noWQ 



1 



Po, 



(15) 



and the temperature is given by 



T = 



P 



k{ne + ni) 
1 

k{ne + Hi) 



(/iAHiunoWo + Pa) - MAmu 



npf \ ' 
ni / 



(^6) 



The position of each zone is found from Eq. (|13|) , using 
the density as the independent variable to solve for a;(ni). 



1 



A[r(n'),n'] 



2A 



B 



dn', 



(17) 



where A is the total cooling rate, and A and B are con- 
stants, determined by the initial values of density and ve- 
locity behind the reverse shock 



(7 + 1) MA \ , .3 

A = — I 1 /iAmu(noUo) , 



2(7-1) V M 



B 



7 



(7-1) 



— - 1 ) (fiArriuVQ + vqTo) 
/i 



(18) 



(19) 



We continue the calculations until the temperature falls 
below 10^ K. Below this temperature the contribution 
to the X-ray emission is negligible. 

As the gas cools, the differential pressure between ad- 
jacent cooling regions causes instabilities, which can create 
secondary shocks within the cooling gas, particularly for 
high shock velocities (Chevalier & Imamura,[1982; Plewa 
& RozyczkaT^S Plewa 1993"1995' Dopita & Sutherland, 
119961 Sutherland et al.,^03). A complete treatment of 
this requires multidimensional hydrodynamical models. 
The specific case discussed in this paper was studied by 
Chevalier & Blondin 1,1995.,) . Although the hydrodynamic 
structure can be severely affected, it is, however, likely 
that the time-averaged spectrum is close to that of a 
steady flow, as was indeed found by Innes H1992|l . As dis- 
cussed by Sutherland & Dopita H199()|l . secondary shocks 
can nevertheless affect the line strengths of the lowest ion- 
ization species. This effect is of minor importance to us, 
since the lowest ionization stages have little or no impact 
on the X-ray emission. Dopita & Sutherland H1996|l also 
conclude that limited spatial resolution in observational 
data has a more serious effect than neglecting thermal 
instabilities. Because several cooling cells and ranges of 
shock parameters are included in one beam, the effect of 
instabilities is reduced. This should certainly hold for the 
cases of interest here, where the size of each cooling cell 
is relatively small. Therefore, we conclude that neglecting 
instabilities in the cooling region is not likely to have se- 
rious implications for the computed spectrum. The main 
effect should be that there should be some spread in the 
reverse shock velocity. The spectrum may then be a super- 
position of several radiative shocks with different veloci- 
ties. This will also be the case for a reverse shock receding 
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into a clumpy ejecta, as seen in e.g. Cas A (Chevalier & 
KirshnerlMi Chevalier & Oishi lMISl Fesen et al.[2005) 
and inferred for some supernovae, like SN 1995N (Fransson 
et al. 120021). 

4. Chemical composition 

The circumstellar shock sweeps up material lost by the 
star before the explosion, which should be close to solar 
composition for Type II supernovae, although CNO burn- 
ing may modify the relative abundances of C, N and O for 
massive progenitors (e.g., Fransson et al. [2005). The re- 
verse shock, on the other hand, propagates into the ejecta, 
sweeping up material, which may have undergone nucle- 
osynthesis. The composition of the ejecta varies with ra- 
dius, because of the successive burning stages in the pro- 
genitor, although hydrodynamical mixing may have mod- 
ified this stratification. The mixing makes it more likely 
that enriched material will be encountered at higher veloc- 
ities. As the reverse shock moves backward into the ejecta, 
it passes through regions of different composition. If the 
progenitor has already lost part or all of its hydrogen en- 
velope, this may occur already during the early phases of 
the supernova. Examples are the Type lib SN 1993J and 
some Type Iln supernovae, like SN 1995N. 

To investigate the effects of composition on the spec- 
trum, we have calculated four models with the same tem- 
perature and maximum ejecta velocity, but different com- 
position characteristics of the burning zones. One is for a 
solar composition, while the other three are characteristic 
for helium, carbon and oxygen dominated regions, respec- 
tively. Abundances of the latter three models were taken 
from the 4H47 model by Nomoto & Hashimoto (1988 ) and 
Shigeyama et al. I|1994|) and are listed in Table [T] The Ni 
abundances of the non-solar compositions, not given by 
the models, were estimated so that the Fe/Ni ratio would 
be the same as for the solar composition. 

5. Ionization structure 

In a medium with negligible external radiation the main 
contribution to the ionization and recombination comes 
from collisions with thermal electrons. Defining the ion- 
ization fraction of ionization stage m of element Z as 
Xrn = n{Zm)/n{Z), we can write the rate of change of 
in a Lagrangian reference frame as 



at 

+0:m+lX,n+l + Cm-lXm-i] , (20) 

where Cm is the collisional ionization rate from ionization 
stage m to TO-l-1, while am is the recombination coefficient 
from ionization stage m to m — 1. Charge transfer, which 
is only important for temperatures below ~ 10^ K, is ig- 
nored for the models in this paper. Charge transfer has, 
however, been included in a later version of the model, but 



was found to have no effect on the spectrum for the cases 
discussed in this paper. The initial state of the gas behind 
the shock front depends on the state of ionization in the 
pre-shock gas, and is discussed in more detail below. We 
follow each gas element as it cools and recombines after 
being shocked. In general, the recombination time will be 
longer than the cooling time, leading to non-equilibrium 
ionization, with highly ionized species being present at 
lower temperatures than in the steady state case. This is 
taken into account by recognizing that each position be- 
hind the shock corresponds to a specific time since the 
gas entered the shock. The time that has passed since a 
particular zone, now at position x, was shocked, is then 
evaluated from t = v{x')~^dx' . 

Eq. (|20|) is solved for the 13 most abundant elements: 
H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni. 
References for the collisional ionization and recombina- 
tion rates are given in the Appendix. 

In our models we neglect all external radiation, assum- 
ing that only shock heating is important for the ioniza- 
tion and emission. Since the radiation from the unshocked 
ejecta has a very low temperature compared to that of the 
shocked gas, this assumption is justified. Because we are 
only interested in the X-ray emitting gas, the effects of 
the X-rays on the cool shell can also be neglected. For the 
optical and UV radiation the reprocessing of the X-rays 
is, however, important l|CF94|) . In addition, the temper- 
ature is low enough, while the density is high enough, 
for ion and electron temperatures behind the shock to be 
nearly in equilibrium (|FL("96|) . We also neglect thermal 
conduction between the zones. This is discussed further in 
Sect.lfOI 

The state of ionization immediately behind the shock 
is set by that of the pre-shock gas. This is in turn de- 
termined by the X-ray photoionization from the shock. 
A self-consistent solution of the two ionization states is 
therefore needed. As usual, the photoionized gas is char- 
acterized by the photoionization parameter C, = L/uer^ 
(e.g.. Tarter et al. ITU^ . Using Eqs. Q and ^ we find 
that 

C ^ 10.4 ( ( \ ' (21) 

Vma -M/ VIOOO km s^V ^ ' 

For a solar composition C, = 12.0(Vtcv/1000 km s~^)^ 
and for an oxygen dominated composition C. ~ 
20.9(T4ov/1000 km s~^)^. The most important point to 
note here is that C only depends on the reverse shock veloc- 
ity and the chemical composition, but not on e.g., the mass 
loss rate parameter, A^. The pre-ionization does there- 
fore not introduce any more parameters into the problem. 
With Vrev ^ 1000 km s'^, implying fcT ~ 1 keV, and 
using the results of Tarter et al. (,1969) for the ionization 
structure based on a free-free spectrum, we expect that 
carbon will mainly be in the form C v-vii, oxygen mainly 
O VI- VIII, and neon mainly Neviii-ix. This is confirmed 
by the models in Fransson H1984|) . 

To calculate the pre-shock state we use the shock spec- 
trum as input to an updated version of the photoionization 
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code in CF94. The result of this is then used as input to 
the shock calculation. If necessary, the new shock spec- 
trum resulting from this can be used for a new calculation 
of the pre-shock structure, etc. In practice, this converges 
after the first iteration. 

6. Emission 

6.1. Continuum emission 

We include free-free and free-bound emission, as well as 
two-photon decay from H-like and He-like ions. The emis- 
sivity (integrated over An) for the free-free and free-bound 
continuum is given by (Sutherland & Dopita. fl993|l 



exffjb = 2.05 X 



(22) 



where the total Gaunt factor Gc is the sum of the Gaunt 
factors for the free-free and free-bound transitions, and is 
written as 



Gc = Gff + Gfh . 



(23) 



Sutherland & Dopita (T^^ and Mewe et al. (1986) 
include a Gaunt factor for the two-photon emission as well. 
For this paper we evaluate the two-photon contribution 
explicitly from the level populations, as described below. 
Any quenching of this at high densities is therefore taken 
into account. The bound- free Gaunt factor is evaluated 
according to the method of Mewe et al. 1)1986(1 . which we 
follow also for the free-free emission. 



Gff = ^ VA2X„(Z) 
n- ^-^ 



(24) 



Z,i 



where Az is the abundance of element Z and z„q is the 
effective charge of the ground state, defined by 

^2 



£^io.r'(^)=^ion(H): 



(25) 



Here, n is the principal quantum number, with the modi- 
fication that we use the calculations of Sutherland H1998|l 
for the individual free- free Gaunt factors, gg. 

For the two-photon emission, we calculate the emissiv- 
ity as 



eA27 



Z,z 



n2^{Z, z)A2^{Z, z)Eo(j){y), 



(26) 



where the sum goes over H-like and He-like ions, z, of all 
elements Z. Further, is the density of the upper level 
of two-photon transitions in ion {Z,z) (either 2s^Si/2 or 
2s^So), Eq is the energy difference between the levels and 
4>{y) is the normalized distribution function of the total 
two-photon emission which comes out at energy y, with 
y = E/Eq. For H-like ions we use the distribution function 
of Goldman & Drake H1981|l . while we use data from Drake 
et al. (|T969 ) for He- like ions. The total emissivity is then 



6.2. Line emission 

In a collisionally ionized, low density plasma the main 
processes responsible for populating the excited levels are 
collisional and radiative excitation and de-excitation. We 
can then write the rate of change of the population rim^k 
of level k of ion m as 



dt 



3<k 

~ ^m,fc GkjTlc ^ ^ Tim.jGjkTic 

]>k j<k 

+ ^ ^ TT-m.jiAjkfijk + GjkTlc) 
]>k 

+ '^■m+lQ^m+l,fe^^o + ^)ri-l?^cC'„i_i , 



(28) 



where Ajk is the transition probability, (3kj is the escape 
probability, am+i,k is the recombination rate from ion 
m -|- 1 to ion m, level k and Cm-i is the collisional ioniza- 
tion rate from ion m — 1 to ion m. Ckj and Cjk are the 
rate coefficients of collisional excitation and de-excitation, 
respectively, with 



kj 



^C,ke' 
9k 



{Ej-Ek)/kT 



3 > k, 



(29) 



where gk and gj are the statistical weights of levels k and 
j, respectively. 

Radiative recombination to individual levels is in- 
cluded for He-like, Li-like and Na-like ions, and dielec- 
tronic recombination to Fexvii and Fexviii. For the ions 
where radiative recombination to individual levels is in- 
cluded we scale the rates to the individual levels with the 
total recombination to the ion, which includes the dielec- 
tronic recombination, thus distributing this contribution 
among the levels. We find that recombination to individ- 
ual levels, where it is included, has little effect on the level 
populations. We assume all collisional ionization to take 
place to and from the ground states. Similarly, recombi- 
nation from ion m to m — 1 is assumed to only occur from 
the ground state, so that am.k = for fc > 0. 

In the Sobolev approximation the escape probability 
is calculated from 



1 



Tkj 



(30) 



where the optical depth of the transition k ^ j, Tkj, is 
given by 



Tkj — flmO'kjjO.s 1 



gkJ^ra.j 
gjl^m^k 



(31) 



Here is the density of ion m, akj the cross section for 
the transition, given by 



(27) <Jkj 



Ajknm.k^Q gj 

S'TTVth gk ' 



(32) 



Nymark et al.: Circumstellar Interaction 



8 

and Cs is the Sobolev length, defined as 
dr 

'Cs = wth3-, (33) 
av 

where Vth is the thermal velocity. The Sobolev approxima- 
tion holds as long as the size of the emitting region is larger 
than the Sobolev length, Cg. In the opposite case we use 
the physical thickness of the region instead of Cg ■ As an 
estimate of dv / dr we take the velocity gradient from the 
shock to where the temperature has fallen to half of its ini- 
tial value. This approximation does not hold for low tem- 
peratures, when the velocity gradient steepens. However, 
the net effect is that we overestimate the absorbing length, 
and thus the optical depth, and since we find that even 
in this case optical depth effects are not important (see 
Sect. l7.4|l . the errors introduced by this method are negli- 
gible. When computing the Sobolev length for lines formed 
in a given shell, we use the thermal velocity in that shell. 

Once the level populations have been found, the vol- 
ume emissivity, ekj, in each line is computed according 
to 

Cfej = nm,kAkjEkj(3kj, (34) 

where Ekj is the energy difference between levels k and j. 

A list of all ions for which level populations are cal- 
culated is given in Table IXTI in the Appendix, as well as 
references for the atomic data. The majority of the atomic 
data have been taken from version 4.2 of the Chianti 
database (Dere et al. nWTI Young et al. l2UU3|l . After this 
paper had been submitted, we tested the code with version 
5.1 of Chianti (Landi et al. I2005|l . but it gave only negli- 
gible effects. To speed up the calculations, only the level 
populations for ions which have an appreciable abundance 
are calculated in each zone. 

The total cooling per unit volume in Eq. (|17|l is calcu- 
lated by adding the contributions from line emission and 
continuum processes 

^ -I- ^ ekj = noniA(ro). (35) 
A k,j 

6.3. Conduction 

As mentioned previously, we have neglected thermal con- 
duction in the cooling region. Previous studies have indi- 
cated that while this is a reasonable approximation for cos- 
mic abundances, a higher metallicity will increase the im- 
portance of conduction (Borkowski et al. 119891 Borkowski 
& ShuUTOT). 

The importance of conduction may be estimated from 
the ratio between the conduction length scale and the cool- 
ing length scale, a ~ icond/-^cooi (Borkowski et al. 11989(1 . 
The conduction length gives an estimate of how far heat 
can be conducted, while the cooling length determines the 
slope of the temperature gradient. Consequently, conduc- 
tion dominates if the conduction length exceeds the cool- 
ing length, i.e. if a ^ 1. We estimate this parameter di- 
rectly behind the shock in the solar and the oxygen domi- 
nated models. For solar composition we find a <C 1, while 



for the oxygen model a ~ 0.02, which justifies neglect of 
conduction even in this case. 

This conclusion is in agreement with Borkowski et al. 
I|1989(l for the solar model, but not for the oxygen model. 
They argue that for a pure oxygen model the cooling rate, 
and therefore also a, is increased by a factor ~ 1000 due 
partly to the higher metallicity, but mainly due to non- 
equilibrium cooling. The latter is, however, an effect of 
the low shock velocities in their models, which means that 
the pre-shock gas is not fully ionized. Therefore, excita- 
tion and ionization of neutrals strongly enhance the post- 
shock cooling, thereby reducing the cooling length. This 
is not the case in our models, where the shock velocity 
is high enough for the pre-shock gas to be fully ionized. 
Instead, non-equilibrium cooling tends to lower the cooling 
rate compared to the equilibrium case, since high ioniza- 
tion stages are present at lower temperatures where they 
cool less efficiently. Thereby the non-equilibrium cooling 
lessens the importance of conduction in our models. We 
therefore feel justified in neglecting thermal conduction in 
the models discussed in this work, but keeping in mind 
that for high metallicities, conduction may become im- 
portant for certain shock parameters, as it would be e.g., 
in supernova remnants where Fe-rich material is encoun- 
tered (e.g., Sorokina et al. l2004|l . It is, however, unhkely to 
change our qualitative conclusions. In addition, any mag- 
netic field present is likely to decrease the conductivity. 

6.4. Line profiles 

The main process responsible for broadening of the emis- 
sion lines is the expansion of the interaction shell, which 
has a velocity of 14 ~ lO'* km s~^. In the thin shell ap- 
proximation we may assume Vs = Vej. For a thin shell at 
radius , the line profile of an optically thin line is flat out 
to the velocity of the shell. For lines which are optically 
thick, the line profile depends on the velocity gradient. 
For V (x r the line is fiat, while for constant velocity it is 
parabolic (e.g., Ambartsumian 1957). A breakdown of the 
Sobolev approximation, as discussed above, may lead to 
double peaked line profiles for the optically thick lines, as 
seen for the Ha line in SN 1998S (Fransson et al.EOISI)- 

In addition, the ejecta will at least during the first few 
years be optically thick to X-rays from the reverse shock. 
Therefore, only the emission from the approaching part of 
the shock will in general be observed. 

Because the line profile in X-rays will usually be dom- 
inated by the instrumental profile, we here use, for sim- 
plicity, a flat, box-like profile in all cases. 

7. Results 

7.1. Model parameters 

The model described here is applicable to supernovae 
where the reverse shock is radiative. Progenitors of Type 
II supernovae have mass loss rates in the range 10^^ — 
lO"'^ M0 yr~^. For a wind velocity = 10 km s^^, this 
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gives A^, « 0.1 — 10, with Type IIPs in the lower range, 
and Type IILs and Type Ilns having A* ^ 1. As we will 
discuss in Sect. 17.31 for r] ^ 10 the reverse shock is then 
radiative to at least a hundred days for solar composition 
and a shock velocity of < 10* km s^^. For steeper den- 
sity gradients, and for helium or oxygen rich ejecta, the 
radiative phase will be longer. The steady state condition 
for the reverse shock is, therefore, at least for the first few 
hundred days, fulfilled for most Type IIL and Type Iln 
supernovae. 

Type IIP supernovae are, because of the moderate 
mass loss rates, in general adiabatic, and are discussed 
in Chevalier et al. H2005(l . Type Ib/c supernovae are 
thought to originate from Wolf-Rayet stars, which have 
much faster stellar winds, with A* ~ 10~^. Therefore, the 
shock is not radiative beyond a few days or even hours, 
and consequently our models are not applicable for Type 
Ib/c supernovae, unless the density gradient is very steep. 
For these supernovae a better assumption is to employ an 
adiabatic model for the temperature and density. Inverse 
Compton scattering of photospheric photons by relativis- 
tic electrons may also contribute to the X-ray emission 
from Type Ib/c supernovae (Bjornsson & Fransson 2003). 
Also inverse Compton scattering by thermal electrons be- 
hind the circumstellar shock may under certain conditions 
contribute fSect. I7.9|l . 

The parameters which characterize a given model are 
the temperature, composition and density immediately 
behind the reverse shock. Further, for a stationary shock 
the density is mainly important for the absolute flux, not 
for the spectral shape, as long as optical depth effects 
and coUisional de-excitation are not important, which is 
in most cases a good approximation (Sect. I7.4|l . As we 
discuss below, departures from ionization equilibrium are 
only important at low temperatures, and do not influence 
the spectrum above 100 eV. Altogether, this implies that 
the shock spectrum is only sensitive to the reverse shock 
velocity and the composition. 

For the models discussed in this paper, we take A, = 1, 
V4 = 1 and s = 2, and a time of 100 days after explosion, 
unless otherwise noted. In order to study the effects of 
different reverse shock temperatures or, alternatively, dif- 
ferent chemical compositions, the density gradient of the 
ejecta is varied according to Eq. Q , as we vary either the 
temperature or the composition. We stress, however, that 
although the density behind the reverse shock varies with 
77, these parameters are only important for the scaling of 
the absolute flux. 

The normalization of the spectrum is given by the fact 
that the total luminosity from a radiative shock is 

Lrov = 47ri?f-pojK-?v 

^ 1.6 X io4i (^-3)(ry --4) ^^ 

vn - 

if s = 2. Because the velocity depends on time as 1/ oc 
^-l/(r,-2)^ this means that the total luminosity of the re- 



verse shock decreases slowly with time, as long as rj is 
constant. However, for the observed flux this is more than 
compensated by the decreasing optical depth in the cool, 
absorbing shell, which lets through more of the emitted 
radiation as the region expands (see Sect. 17. 8|) . Note also 
that approximately half of the above luminosity is ab- 
sorbed by the ejecta (Sect. EHI IZ3- 

7.2. Shock structure 

Radiative shocks occur in a wide range of circumstances 
in the interstellar medium, and have been discussed by a 
number of authors. The first numerical models of cooling 
shocks were made by Cox (QW2|, who considered a pla- 
nar shock with a velocity of 100 km s~^ and applied the 
results to the Cygnus Loop. Later, improved models were 
developed by Raymond (jMHI), Dopita ifWHI IHT?! and 
others (see also Dopita & Sutherland 1 1 9961 and references 
therein) . These models were intended for such diverse phe- 
nomena as old supernova remnants and narrow line emis- 
sion regions in active galactic nuclei, but are in general 
not applicable to young supernovae. The reason for this is 
that for the low densities in supernova remnants the shock 
is not radiative above ~ 200 km s~^. 

The effects of composition on the properties of the 
cooling gas in radiative shocks have been studied by Itoh 
|1981a |0-981b, 1986, 1988), Borkowski et al. ijTOg^ and 
Borkowski & ShuU P99b|l . The main difference between 
these studies and the present one is again the lower ve- 
locity, 5s 200 km s~^, and density in their models, which 
is the reason for the somewhat differing results, e.g. as 
regards the role of conduction, as mentioned earlier. 

The hydrodynamic evolution of fast shocks evolving in 
a dense medium, applicable to AGNs and young supernova 
remnants, was calculated by Plewa & Rozyczka (|1992f) and 
Plewa 11993, 1995). While these models could indeed be 
applicable to supernovae if the composition is close to so- 
lar, their treatment of the cooling in a time-independent 
formalism makes it less applicable to the reverse shock. As 
we discuss below, a higher metallicity, as appropriate for 
the reverse shock, leads to a stronger cooling, and thereby 
a need for a time-dependent calculation of the ionization 
balance and cooling. They also treat the emission in a 
simplified manner, assuming coUisional ionization balance 
and computing the emission in the cooling region from 
given line intensities. In contrast, we compute the emis- 
sion from transitions in multi-level ions, and couple the 
resulting cooling to the hydrodynamics resulting in de- 
tailed and self-consistent spectra. 

An important feature of high velocity shocks is the 
creation of ionizing photons, which ionize the pre-shock 
gas and the cooling, shocked gas (Sect. EJ. If neutrals 
are present, some time is needed for the shocked gas to 
become ionized. As the shocked gas approaches ioniza- 
tion equilibrium corresponding to the post-shock temper- 
ature, the non-equilibrium ionization state creates strong 
line emission, which causes rapid cooling. For very fast 
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shocks {Vs ^ 500 km s like the shocks considered in 
the present paper, the gas is, however, nearly fully pre- 
ionized (Sect.EI Fransson [ro51lCF94(l . 

When ionization equilibrium has been reached, the gas 
cools slowly, with the ionization state adjusting to equi- 
librium at the local temperature more rapidly than the 
cooling occurs. This is where most of the emission occurs, 
because of the high temperature. When the local temper- 
ature falls below that corresponding to the maximum of 
the cooling function, around 10^ K, non-equilibrium ion- 
ization sets in. The gas cools rapidly, until the temperature 
falls below 10^ K, where photoionization heating balances 
cooling. 

The gas close to the contact discontinuity has a tem- 
perature 10^ K and a density of - (lO'^ K/IO^ K) - 10^ 
times that of the post-shock gas, i.e., 10^° — 10^^ cm^^^, 
and absorbs a large fraction of the X-rays. fSect. I7T5|) . A 
large fraction of the optical and UV emission originates in 
the cold, dense shell l|CF94|l . 

For a given reverse shock velocity and composition, the 
post-shock temperature will be lower the steeper the ejecta 
density gradient is, but is typically ^ 10^ K. This is close 
to the minimum of the cooling function, so in models with 
a temperature lower than ~ 10^ K a cooling instability will 
set in almost immediately, whereas models with a higher 
temperature will take longer before cooling starts, i.e., the 
interaction region will be thicker. 
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Fig. 2. Density (dashed lines) and temperature (solid 
lines) behind the reverse shock as a function of distance 
from the contact discontinuity in a model with V4 = 1, 
s = 2, A, = 1, and To = 1.0 keV at 100 days after the 
explosion and for three compositions. In these models the 
radius of the interaction region is Rg = 8.64 x 10^^ cm. 
The vertical lines indicate the position of the shock with 
respect to the cool shell in the three cases. 



the distance from the end of the cooling instability, rather 
than from the shock. 

An analytical estimate of the thickness of the cool- 
ing region may be obtained from the cooling time and 
the velocity of the post-shock gas, Kov/4, so that D — 
icooiKcv/4 where 



^cool — 



3kT{nc 



3kT 



2Aneni 



2(1 - fi/fiA)niA 



(37) 



The equilibrium cooling function may be fit by a func- 
tion of the form 



A = AT, 



6 



6 ' 



(38) 



where Tq — K and the fitting parameters A, B, 

a and /3 depend on the composition of the gas, and are 
given in Table [3 for the compositions considered here. In 
the range 10'°^ — 10® K, which are the temperatures con- 
tributing to the X-ray emission, the fits are in general valid 
to within a factor of 2. The exceptions are the carbon and 
oxygen models, for which the cooling below 10^ K varies 
strongly. The fits to the cooling functions are plotted in 
Fig. 13 together with the cooling functions from our model 
calculations, assuming ionization equilibrium. 
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Fig. 3. Cooling curves for four compositions. Cooling 
functions from the model calculations for each composi- 
tion, assuming ionization equilibrium, are plotted with the 
corresponding analytical fit, with the fitting parameters 
given in Table 12 

For the temperatures of interest, only the first term of 
the cooling function in Eq. H38I) contributes. Combining 
this with Eq. (|S7|l . and using Eqs. and © for the 
temperature, we can express the thickness of the region as 



In Fig. [3 we show the density and temperature struc- 
ture for our standard model with V4 = 1, A, = 1 and 
To = 1.0 keV at 100 days after the explosion, and for three 
compositions. To show the drop in temperature and ion- 
ization close to the cool shell more clearly, we plot this as 



D 



1.45 X 10-** AiAM^+" 
A {fiA - 
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(7?-3)(77-4)(7y-2) 



3+2a 



5+2a 



(39) 



where the time is in days since the explosion, and we have 
set s = 2. The composition dependent constant Dq is given 
in Table 12 for the four compositions considered here. 

For the model parameters above (V4 = 1, A, = 1, Tq = 
1 kcV), the analytical estimate gives D ^ 2.1 x 10^'' cm, for 
solar composition, while the model calculations give D ~ 
5.3 X 10^^ cm. The factor of four difference results from 
the approximate form of the cooling function in Eq. (|38|l 
and the fact that the cooling increases as the temperature 
drops, leading to a collapse of the region. For the non-solar 
models the analytical estimates give D ~ 3.5 x 10^'^ cm 
for the helium model, D ^ 2.8 x 10""^^ cm for the carbon 
model and D ^ 6.7 x 10^^ cm for the oxygen model. The 
numerical results for these models are D ^ 7.8 x 10^^ cm, 
D ~ 9.4 x 10^^ cm and D 3.1 x 10^^ cm, respectively. In 
all cases the thickness of the cooling gas is much smaller 
than the shock radius, justifying the thin shell and sta- 
tionary approximations. 

The ionization structure of Fc in the cooling region 
behind the reverse shock for the solar composition model 
is plotted in Fig. 2| We see that Fe is never fully ion- 
ized, because of the comparatively low shock temperature. 
Immediately behind the shock front Fe xxi is most abun- 
dant, with noticeable fractions of Fexvii-xxv. These ions 
dominate the bulk of the cooling gas, but further down- 
stream also lower ionization stages become sufficiently 
abundant to affect the emission. 
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expansion time. In this case we can calculate the tempo- 
ral evolution assuming steady state with given values of 
no and Tq for each epoch. For the hydrodynamic equa- 
tions the relevant timescale is the cooling timescale, given 
by Eq. |T7|l. 

We check the assumption of a stationary flow by esti- 
mating tcooi in all zones. Using Eq. (jSJ for the ion density 
and s=2, we can express the ratio of the cooling and ex- 
pansion timescales as 



cool 



6.7 X IQ-^^nlfiV^td 



t (MA-M)(^-2P(??-3)(77-4)A(r)A, 



(40) 



where is the time in days since the explosion. The cool- 
ing time is thus determined by the shock velocity, Vs , and 
the ratio t^/ A^. 

In Fig. we plot this relation for s — 2, three differ- 
ent values of 77 and three of the chemical compositions in 
Table ^ For a particular value of A^, the reverse shock is 
radiative at a given time, id, if the velocity of the outer 
shock is below the curve corresponding to the ejecta den- 
sity profile. We see that the higher the value of rj, the 
larger the range of shock velocities for which the reverse 
shock becomes radiative early on. The same is true as we 
go from solar composition to more metal rich gas. 
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Fig. 5. Ratio of cooling time (icooi) and age (id) of the 

shocked ejecta. The lines represent the limiting velocities 

for which the cooling time behind the reverse shock is 

equal to the expansion time. If the velocity of the circum- 
Fig. 4. Ionization structure of Fe behind the reverse shock ^^^u^^ ^^^^^ ^ ^.^^^ 1^^^^ ^^^^ ^^^^^^ ^l^^g 

for the solar composition model in FigM The vertical hue ^^^^^^ corresponding to the relevant density slope, ry, 

indicates the position of the shock. ^^^^^^^ ^^^^^ radiative. If it is higher, the shock is 

adiabatic. 



7.3. Cooling and ionization timescales 

The main assumption of our model is that the flow be- 
hind the reverse shock is stationary, i.e., that the time 
needed for a process to occur is short compared with the 



Due to the nearly constant pressure behind the reverse 
shock, ric (X in the post-shock region. The density in 
a model with To ~ 10^ K will therefore increase by a 
factor ~ 10^ before reaching To ^ 10^ K. Therefore, if the 
assumption of stationarity holds immediately behind the 
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reverse shock, it will be a good approximation in all of the 
post-shock gas. 

Although it has been found that in a supernova rem- 
nant the ionization structure must be modeled by a 
time-dependent formahsm (e.g., Kaastra & .Tansen IT99HI 
Sorokina et al. 2004), it is not clear if this is also the 
case for a young supernova, where the densities are much 
higher. We will therefore next explore the differences be- 
tween steady state and time dependent calculations for the 
ionization balance for a variety of parameters and compo- 
sitions (Table [U. 

For the ionization balance, the requirement is that 
the collisional ionization parameter t^ei = ngtrei be suf- 
ficiently small compared to the expansion and cooling 
timescales. The relaxation time for a particular ion, treh 
is the shortest of the ionization and recombination times 
(e.g., Mewe HM^ . given by U^n = [C{Tc)nc]^^ and 
tree = [a(7c)no]~^, where C{Te) and a{Te) are the col- 
lisional ionization and recombination coefficients, respec- 
tively, for a given ion at a specific temperature. There 
is, therefore, a specific, temperature dependent relax- 
ation parameter for each ion, which can be expressed 
as troi = min[l/C(rc), 1/q:(T'c)]. In Fig. |H| we show how 
this parameter for Fe ii-xxvi varies with temperature for 
T ^ 10^ K. At high temperatures it is given by the ion- 
ization timescale, while at lower temperatures it is deter- 
mined by the recombination time. The relaxation param- 
eter depends only on atomic physics. 

Let us first consider the ratio of troi to the expansion 
timescale, t. For a given temperature, a limiting density 
above which steady state applies may be found from 
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net > irol(T) 



(41) 



Using the ion density from Eq. Q, writing the electron 
density as rto — '^ion(AiA — ^-i^d assuming s=2, this 

relation translates into the requirement that the ejecta 
density gradient satisfy the following relation 



(7/-3)(7/-4) ;> 0.12 



10"^ days cm.~^ J \100 days 



(42) 



The highest values of the relaxation parameter in Fig. 
are those of Fexxv and Fexxvi at T = 10^ K, where 
frei ~ 10^ days cm~'^. This is also among the highest val- 
ues of iroi for any ion, and can therefore be used as an 
upper limit to troi. For the set of parameters which we use 
in this paper (A* = ¥4 — 1, see Sect. l7.l|l we find that this 
relation is fulfilled for both the solar and higher metallic- 
ity models up to several hundred days for all interesting 
values of rj. Only for shocks with very high velocity and 
low mass loss rate is the expansion timescale shorter than 
the ionization timescale. 

For steady state ionization to hold, the relaxation time 
for ionization must also be shorter than the local cooling 
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Fig. 6. The lines represent the collisional ionization pa- 
rameter, tici = trc\nc for Fc, with t^^i h^ days. 

time, frci < icooi- The cooling time decreases rapidly as 
the temperature drops in the post-shock gas, and a shock 
for which steady state holds at the shock front may be 
out of equilibrium further downstream. In terms of the 
collisional ionization parameter and the cooling function 
we get 



(43) 



Using the first term of the cooling function in Eq. (|38|) . 
for Te ^ 10^ K, we can express the condition for steady 
state ionization as 



trcl < 



4.8 X 10" 



HA A 



Tgi+" days cm" 



(44) 



This requirement should be checked for each ion at the 
temperatures where the ion is abundant. Note that this 
only depends on the temperature and composition, but 
not on the density. 

As an example we consider Fe in the solar composi- 
tion model, for which the collisional ionization parameter 
is shown in Fig.El Using trci ~ lO'' days cm^"^ as an upper 
limit to trci, as discussed above, we can determine a lower 
limit to the temperature, ~ 5 x 10^ K, above which steady 
state can be assumed to hold. Below this temperature time 
dependent effects are expected to become important, al- 
though the fact that t^^i varies with the temperature, and 
that the value used here is an upper limit, means that 
steady state will be a reasonable approximation also at 
somewhat lower temperatures. The corresponding require- 
ment for the helium model is T ^ 1.8 x 10^ K, while for 
the carbon and oxygen dominated models T ^ 1.1 x 10^ K 
and T ^ 3 X 10^ K, respectively. 

From this we expect steady state to be a reasonable 
approximation for the X-ray emitting part of a solar com- 
position and a helium dominated shock. The carbon and 
oxygen dominated shocks, on the other hand, are expected 
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to be out of equilibrium already at the shock front. For 
these shocks a fully time dependent ionization balance cal- 
culation is required, as we indeed do in this paper. 

To illustrate this quantitatively we show in Figs.[7|and 
IHl the ionization structure of Fe both for steady state and 
time-dependent ionization in our standard model with so- 
lar composition and oxygen dominated composition, re- 
spectively. We see that in the solar case the state of ion- 
ization behind the shock follows steady state until the gas 
has cooled to a few times 10^ K. After that, recombination 
begins to lag behind coohng, the effects being noticeable in 
Fevi and lower ions, which is what we would expect from 
the estimate above. Since this occurs at such low temper- 
atures, it has no noticeable effect on the X-ray emission. 

For the carbon and oxygen dominated models non- 
equilibrium effects are important already for Fexxi at 
7 X 10^ K, as we expected. As extreme examples, we note 
in Fig. IHl that the concentrations of FeXVll and FeXVlll 
stay high at considerably lower temperatures than in the 
steady-state case, due to the long recombination time com- 
pared to the cooling time. Except for A ^ 100 A, the effects 
on the spectra are, however, minimal. 

From this discussion we conclude that the equilib- 
rium assumption is reasonable, except for low temperature 
shocks, and shocked gas highly enriched in heavy elements. 
This is useful since it allows a much simplified calculation 
of the shock emission, without a fully time dependent for- 
malism. 




10= 10^5 10^ 

Temperature (K) 



Fig. 7. Ionization structure of Fe behind the reverse shock 
as function of temperature. The thin lines are the steady 
state model, and the thick lines are for a time-dependent 
ionization. The input parameters are the same as for the 
model in Figs. 13 and 0] and the composition is solar. 



7.4. Optical depth and density effects 

Radiative transfer is taken into account by means of the 
escape probability formalism, described in Sect. but 




Temperature (K) 

Fig. 8. The same as in Fig.|3 but for an oxygen dominated 
composition. 

turns out to have little effect. To understand the reason 
for this, we estimate the optical depth, r, for typical pa- 
rameters. Combining Eqs. H31() and H32() . and using the 
total size of the region, D, instead of the Sobolcv length, 
as is more realistic, 

T ^ 4.9 X lO-i2^ei?A(Z)X„-4^, (45) 

where A{Z) is the abundance of element Z , Xm is the ion- 
ization fraction of ion m, and the energy of the transition, 
Ekj, is in eV. 

As particularly interesting examples we consider two 
of the strongest lines, the He-like Ovii A21.60 and the 
H-like Oviii A18.97 resonance lines, with / = 0.69 and 
/ = 0.27, respectively. For collisional ionization equilib- 
rium the maximum abundance of O Vii occurs at T '-^ 
8 X 10^ K, where Xm 1, and for O viii at T ~ 2.5 x 10*^ 
K, where Xm ^ 0.5. If we now write 

T = ToA{Z)N2u (46) 

where N21 is the total column density in units of 
IG^i cm^^ we find that tq = 6.6 x 10^ for Ovii A21.60, 
and To = 6.4 x 10^ for Oviii A18.97. 

The column density of the cooling gas is 

and therefore depends only on the composition and reverse 
shock velocity, but not on the mass loss rate or shock 
radius. For solar composition and Tq ^ 10^ K we get N ^ 
10^^ cm~^. Further, A{0) — 6.9 x 10^"', and we therefore 
obtain r 4.4 for Ovii A21.60, and r = 0.44 for Oviii 
A18.97. 

For an oxygen dominated composition we get in the 
same way t 5.6x 10"^Af2i for O vii A21.60. Because of the 
stronger cooling the radiating region is also much thinner 
(Sect. I7.2|l . Our model calculation yields a thickness of 
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4 X 10^"'^ cm and a column density N21 ~ 6.3 x 10"'^ for 
the same input parameters as above. Therefore, the higher 
oxygen abundance results in an optical depth of O Vii that 
is comparable to that of the solar model. Some of the 
strongest lines may therefore be marginally optically thick. 
The great majority have, however, t ^ I. 

A consequence of an optical depth higher than unity 
is that coUisional de-excitation or transfer between dif- 
ferent L and S states may become important. The ratio 
of collisional de-excitation to radiative escape is given by 
ricCkj / f}k]Akj ~ UcCkjTkj/Akj, for Tkj ^ 1. 

For the O vii resonance line the most important de- 
struction channels are either by collisional de-excitation, 
or by a collisional transition to the 2s ^Sq level, followed 
by a two-photon transition. At 10® K the collisional de- 
excitation rate is ap^-^is = 2.4 x 10~^° cm'^ s~^, 
while the radiative transition rate is ^2s3Pi^isiSo ~ 
3.2 X 10^^ s^^. Therefore, there is a critical density, 
ncrit,2s3Pi_isiSo = 2.2 X 10^^N^^^A{O)-'^ cm'^, above 
which collisional destruction dominates radiative de- 
excitation. The forbidden 2s 6*1 Is ^6*0 line has 
C2s35i^isi5o = 4.1 X 10"" cm^ s"i and vla^ssi^is iSo = 

1.1 X 10^ s-\ and n„it,2s3Si-isiSo = 2.7 x lO^^ cm'^. 
Transitions to the singlets occur with C'2s!>Si—f2s'^Sa ~ 

1.2 X 10~^° cm"^ s~^, so conversion to two-photon emission 
occurs at a somewhat lower density, ~ 10^"^ cm""^. To sum- 
marize, conversion to two-photon emission or destruction 
by collisional de-excitation will only be important above 
~ lO^'^ cm"'^, independent of the optical depth. 

Before this occurs, collisional excitation from 2s^Si 
to 2p^P2 wiU for n ^ ncrit,2s3Si-2p3P2 = 7 x 10^" cm"^ 
convert radiative decays in the A22.10 2s^So Is^Sq 
line to the A21.80 2p Is^Sq line. Because of limited 

instrumental resolution, as well as velocity smearing, the 
individual lines will, however, be difficult to resolve. The 
total flux in the A21.60 — 22.10 lines will therefore not be 
affected below ^ 10^^ cm~^. 

To check the dependence on the density and optical 
depth we have made calculations for two sets of parame- 
ters, one set with A* = 1 at 100 days, and one set with 
^4* = 10 at 10 days, and both with 77 = 10, s = 2 and 
Tq — 1 keV. The density behind the shock in the lat- 
ter calculations is therefore a factor of 10"^ higher than in 
the 'standard' calculations below. This range covers most 
cases of interest. 

After adjusting for the factor of ten difference in flux, 
because of the difference in mass loss rates, we find that 
only the He-like Cv A41.47 2^ So I'^Sq and A40.73 
2^P2 l^So lines, and the corresponding Ovii A22.10 
and A21.80 lines are affected by more than 50%, as pre- 
dicted from the simple arguments above. Other lines are 
largely unaffected by the difference in density. 

We therefore expect that neither density nor optical 
depths should affect the fluxes appreciably, except for 
some He-like ions. Even in these cases the effect is likely 
to be hidden by the instrumental resolution or velocity 
broadening. 



7.5. Variation with reverse shock velocity and 
composition 

From the discussion in Sect. 17.31 we find that as long as 
the steady state requirement is fulfilled, the density behind 
the reverse shock is of secondary importance for the spec- 
trum, except as a normalization parameter fEa. I36|l . We 
therefore now discuss the effects of varying the remaining 
parameters, the shock temperature and the composition. 

7.5.1. Reverse shock velocity 

The temperature of the post-shock gas depends on the 
velocity of the reverse shock. In Fig. IHl we compare the 
spectra from four models with solar composition, where 
the reverse shock temperature, Tq, is 0.3 keV, 1.0 keV, 3.0 
keV and 10.0 keV, respectively. For a given shock velocity, 
each temperature corresponds to a particular value of the 
ejecta density gradient given by Eqs. Q and jnj. Here we 
have fixed the maximum velocity of the ejecta at V4 = 1. 
The above reverse shock temperatures then correspond to 
density gradients 77 = 22.0 for Tq — 0.3 keV, 77 = 8.3 for 
To = 3.0 keV and 77 = 5.5 for To = 10.0 kcV. The ejecta 
velocity is only important for the velocity broadening of 
the lines. To show the distribution of the flux with wave- 
length more clearly, we plot the luminosity per logarithmic 
wavelength interval, ATa- 

The most obvious difference between the four models is 
the high energy exponential cutoff. Further, a higher tem- 
perature leads to a higher degree of ionization immediately 
behind the shock. However, because the gas is cooling, low 
and intermediate ionization stages are present even in the 
hottest model. In the models with To <; 3 keV, the hot gas 
behind the shock produces little or no line emission, while 
contributing strongly to the continuum emission, which 
overwhelms the line emission from the cooler region. 

Fig.llOlshows the ratio of the strongest lines to the fiux 
of the strong Oviii A18.97 line as a function of Tq for so- 
lar composition. The left hand panel shows the strongest 
lines of the metals up to Ca, while the right hand panel 
shows the strongest Fe lines. Note that these relative lu- 
minosities do not include the continuum contribution. It is 
seen from this that at low temperatures the line emission 
is dominated by C, N and O, while the H- and He- like ions 
of Mg, Si and S become stronger for To ^ 1 keV. Lines 
of especially Fexix, Fexxi, Fexxiii, Fexxiv and Fexxv 
dominate successively the cooling for Tq 0.4 keV. 

The most conspicuous features in the To = 1 keV and 
To = 3 keV models are the Fe complexes between 10-15 A 
and around 100 A. The complex at 10-15 A consists of 
emission from Fexxi-xxiv, mixed with Oviii and Neix- 
X, while the 100 A feature is a blend of lines from lower 
ionization stages of Fe, Fe xviii-xxii, as weU as O vii-viii, 
Neviii and Nixxiii-xxv. These features are also present 
in the To — 10 keV model, but have a smaller equivalent 
width due to the strong continuum. In the Tq = 0.3 keV 
model the Fe features are practically absent, because the 
high ionization stages of Fe are no longer present. This can 




Fig. 9. Spectra with different reverse shock temperatures, To = 0.3, 1.0, 3.0 and 10.0 keV, ali for solar composition. 



be seen in Fig.|3 which shows that at a temperature of ~ 
3.5 X 10^ K Fe XVIII is the most abundant ion, while higher 
ionization stages of Fe are practically non-existent. Lower 
ionization stages of Fe have few emission lines at these 
energies, and the line emission in this model is therefore 
taken over by lighter metals. 

Another line complex is seen at 175-200 A, which con- 
sists mainly of lines from Fex-xi, with contributions from 
intermediate stages of O, Ne, Ar, Ca and Ni. This feature 
is strong in the Tq = 0.3 keV model, but weakens with 
increasing shock temperature. 

The emission at A ^ 10 A comes from H-like and Hc- 
like stages of Ne, Mg, Si, S, Ar, Ca, Fe and Ni. From Fig.lTUI 
it is apparent that these lines become increasingly stronger 
as the temperature increases. In the Tq — 3 keV spectrum 
the Fe K emission is the strongest, with the Fexxv 1.85 A 
line dominating, as well as weaker Fexxv AA 1.46, 1.57- 
1.58 lines. In the Tq = 10 keV model the Fe K hues at 
1.4-1.6 A are strong, and now include Fexxv-xxvi, as 
well as a contribution from Nixxvii at 1.60 A. 

In the To — 0.3 keV model we also note the O viii 
absorption edge at 17 A, which is not seen in the other 
models. 



7.5.2. Effects of composition 

As described in Sect. ^ we have calculated four mod- 
els with different chemical composition, one with solar 
composition, and three corresponding to the most impor- 
tant burning zones in the progenitor (Table In Fig. 1111 
the spectra of the four models are shown. In all cases 
To = 1 keV. As we explained in Sect. 17.11 we assume 
that V4 = 1 in all cases, and to get the same Tq for all 
compositions we have varied rj. For V4 = 1 a reverse shock 
temperature Tq = 1.0 keV corresponds to 77 = 13.0 for the 
solar model, rj = 18.3 for the helium dominated model, 
f] = 20.2 for the carbon dominated model and t] ~ 20.7 
for the oxygen dominated model. 

The most conspicuous features in the helium and solar 
models are the Fe line complexes around 10 A and 100 A. 
These features are enhanced by lines from Neix-X near 
10 A, and Ne VII- VIII and Ovii near 100 A, otherwise Fe 
dominates these features. The higher abundance of C, O 
and Ne in the helium model means that these elements 
dominate the cooling, which can also be seen from the 
strength of the lines, particularly from Cv-vi and Neix- 
X. Therefore, lines from heavier metals are stronger in the 
solar model than in the helium model, in spite of their 
lower abundance. This is especially apparent below 10 A, 
where the line emission in both models comes from H-like 
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Fig. 10. Luminosities of the strongest lines relative to O viii A18.97 as function of the temperature of the reverse 
shock. The flux weighted wavelengths of the different Hues labeled are Cv A41.5, Cvi A33.7, O vii A22.2, Neix A13.7, 
Nex A12.1, Mgxi A9.3, Mgxii A8.4, Sixiii A6.7, Sixiv A6.2, Sxv A5.1, Sxvi A4.7, Fexvii A17.1, Fexviii AA14,6, 
17.7, Fexix AA14.0, 16.8, Fexx A15.5, Fexxi A14.2, Fexxii A13.2, Fexxiii A12.2, Fexxiv AA8.4, 11.4, Fexxv A1.87, 
Fexxvi A1.78. 



and He-like Mg, Si and S and He-like Ar, all of which are 
stronger in the solar model than in the helium model. 

In the carbon and oxygen dominated models the Fe 
features are more or less absent. This is because of the 
strong overabundance of O and Mg, as well as Si and S, 
which take over the cooling in the high energy part of the 
spectrum. Therefore, although the Fe abundance is higher 
than in the other models, the Fe lines are overwhelmed by 
the stronger emission from lighter elements. Consequently, 
instead of the Fe features near 10 A and 100 A, the oxygen 
model has strong lines from a number of ions distributed 
over the whole spectrum, while the carbon model has few, 
but strong lines. 

The line emission of the carbon model is dominated 
by lines from Cv-Vl, Ovil-Vlll and NeVll-X. The oxygen 
model, on the other hand has only weak C emission and no 
Ne lines at all as a result of the low C and Ne abundances. 
Instead, it is dominated by emission from O VII-VIII, mixed 
with strong lines from H-like and He- like Mg, Si and S, and 
lower ionization stages of the same elements between 10 - 
100 A. 



The continuum of the carbon and oxygen models is 
considerably steeper at long wavelengths compared to 
the solar and helium models. This reflects the increasing 
importance of bound-free and two-photon contributions 
compared to the free-free emission in the metal rich mod- 
els. As a result both the carbon and oxygen models show 
strong absorption edges, which are not apparent in the so- 
lar model. The most prominent is the O IX edge at 13.8 A, 
which is also present, albeit weak, in the helium model. In 
addition, we see strong edges from Nexi at 51.7 A, and 
C VII at 24.8 A in the carbon model, and from Mgxiii at 
7.3 A and Mgxiv at 6.5 A in the oxygen model. Also the 
absorption edges of Neviii and Si VII, which coincide at 
62.0 A, are apparent in the oxygen spectrum. In both the 
carbon and oxygen models the signature of two-photon 
emission from O Vli and O VIII is apparent, most clearly 
in the carbon dominated model, where it is not obscured 
by line emission. 
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Fig. 11. X-ray spectra produced by four models with different composition. All models have V4 — 1 and Tq = 1.0 keV. 
The chemical abundances in each zone are listed in Table [H 



7.6. Effects of pre-ionization 

As mentioned previously the state of ionization immedi- 
ately behind the shock is determined by that of the pre- 
shock ejecta gas. For the 1 keV case most pre-shock oxygen 
is in Ov-vii, while iron is found as Feviii-xi, both for so- 
lar and oxygen dominated composition. Only after ~ t^-d 
does the ionization reach the state corresponding to the 
local temperature. 

In Fig. El we show the difference between the emitted 
spectrum when pre-ionization is included, and when we 
assume that ionization equilibrium is reached instantly, 
for a solar composition and for an oxygen dominated gas. 
In the solar mixture the difference is barely noticeable, ex- 
cept for very short wavelengths, and even there it is neg- 
ligible. In the oxygen dominated gas, on the other hand, 
the difference is pronounced for wavelengths shorter than 
~ 15 A. The reason for this is the much shorter cooling 
time in the oxygen dominated model compared to the so- 
lar model. The steepening continuum emission at short 
wavelengths in the oxygen dominated model is caused 
by the dominance of free-bound and two-photon emis- 
sion over free-free emission, since both free-bound and 
two-photon processes are strongly affected by the lower 
ionization state in the post-shock gas when pre-ionization 



is included. For the line emission, the hydrogen-like lines 
are most affected by pre-ionization, as expected. We note 
e.g., that the hydrogen-like Mgxii line at 8.4 A is a factor 
~ 3 lower when pre-ionization is included. The helium-like 
Mgxi line at 9.3 A, however, is hardly affected at all. This 
is as expected, since the effect of pre-ionization is to re- 
duce the abundance of high ionization stages close to the 
shock front. 

7.7. Comparison with single-temperature models 

When X-ray observations of supernovae are analyzed, the 
emitting region is usually assumed to have a single tem- 
perature. In some cases two temperatures are used, to ac- 
count for the different temperatures of the two shocks. As 
an example, Uno et al. H2UU2|I fit the ASCA observations 
of SN 1993J with two components, with different temper- 
atures and absorption column densities. However, while 
a single temperature fit is a good approximation for the 
outer shock, which is adiabatic, it is far from sufficient for 
reproducing the emission from a radiative shock, as we 
will show below. 

Another interesting example is the Chandra obser- 
vations of SN 1993J after 8 years, which could only be 
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Fig. 12. Comparison between spectra where pre- ionization is included (thick hne) and spectra where ionization cqui- 
hbrium is assumed to be reached instantly in the post-shock gas (thin line) . The left hand panel is for solar composition 
and the right hand panel is for an oxygen dominated composition. In all cases V4 = 1 and Tq = 1.0 keV. 



fit with two low temperature components (0.35 keV and 
1.01 keV) and one high temperature component (6.0 keV) 
(Swartz et aL l^UnSjl . Similarly Schlegel et al. if^HMjl find, 
when trying to fit spectral models to recent Chandra 
observations of SN 1978K, that while a two component 
model gives the best fit of the models tested, the fit even 
for this model was not satisfactory. This indicates that 
more detailed models are necessary in order to fit the ob- 
served spectra. In particular, it is important to realize that 
the relative contribution of the different temperature com- 
ponents can not be treated as free parameters, but are 
fixed by the hydrodynamic structure of the shock as it 
cools. 

In Fig. El a spectrum produced by our combined 
hydro-/spectral code, with Tq = 1.0 keV, is compared to a 
single temperature spectrum with a temperature equal to 
the post-shock temperature in the whole region between 
the reverse shock and contact discontinuity. The spectra 
are normalized to have the same total luminosity between 
1-10 keV. Some of the most prominent lines in each 
model are given in Table|31 where the fluxes relative to one 
of the most prominent lines in the spectra, O viii A 18.97, 
are listed for these two models, as well as for a single tem- 
perature model with T = 0.7 keV. 

As expected, the continuum emission in the cool- 
ing model is higher at low energies, and there are more 
lines from lower ionization stages than in the single- 
temperature model. In particular, we see that the line 
complexes around 10 A and 100 A are broader in the full 
model than the single-temperature model, a result of the 
wider range of ionization stages contributing to these fea- 
tures. At energies below 10 A, the line emission in both 
models is dominated by H-like and He-like Mg, Si and S, 
the H-like lines being stronger in the single-temperature 
model than in the full model. At wavelengths longer than 
~ 15 A the line emission is stronger in the full model 
than in the single-temperature model, due to the presence 



of colder gas. The lines in the single-temperature model 
mainly come from Fexx-xxiv, together with weaker lines 
from Cvi, Oviii and Nex, while the high temperature 
does not allow lower ionization stages. In the cooling 
model the same lines are present, but also lines of lower 
ionization stages like C v, N vii and O vii, as well as Si IX- 
XIII, Sxi-xiv, Arxiv-xvi, Caxv-xviii and Fexv-xvi. 

Quantitatively, the difference can be seen from Table|51 
The line ratios of, e.g., the H-like lines of the different ele- 
ments differ by large factors between the two models. In re- 
ality one would not use the temperature behind the shock 
for a single-temperature model, but instead try to find 
the temperature which gives the best fit to the observed 
spectrum. Because of the contributions from cooler zones, 
the best-fit spectrum will correspond to a lower temper- 
ature than the shock temperature. When attempting to 
fit the shock model with Tq = 1.0 keV shown in Fig. [T^ 
with a single-temperature model, we find the best overall 
fit to be a single-temperature model with T = 0.7 keV. 
Although there is a slight improvement compared to the 
single-temperature T = 1 keV model, the line ratios differ 
by at least a factor of two from the shock model. In ad- 
dition, the 'best-fit' temperature is too low. To reproduce 
the line emission from both high and low ionization stages 
a combination of cold and hot gas is therefore necessary. 
Note, however, that a multi-component model is not suf- 
ficient, because the contributions from each of these com- 
ponents have to be calculated self- consistently, as in the 
shock model. The large differences between these models 
are of obvious importance when determining abundances, 
and a too simplified analysis can lead to large errors in 
the results. 



7.8. Absorption by the cool shell 

When the reverse shock is radiative, most of the emission 
from the reverse shock is absorbed by the cool shell close to 
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Fig. 13. Comparison between a single temperature spec- 
trum with T — 1.0 keV (thin Hne) and a cooUng shock 
spectrum created by our model, with Tq = 1.0 keV (thick 
line), both with solar composition. The spectra are nor- 
malized to have the same total flux in the range 1-10 keV. 

the contact discontinuity. At early epochs, only the most 
energetic emission passes through, while lower energies are 
completely absorbed. As the emitting region expands, the 
column density in the cool shell drops, and the shell be- 
comes transparent, first at high energies, later at lower 
energies. The absorption can be taken into account by 
calculating the emerging luminosity at wavelength A as 

Lx - LA,cme-"(^) , (48) 

where the subscript 'em' refers to the emitted luminosity, 
T — A'coio'(A), and (t(A) is the absorption cross section at 
wavelength A, given by 

(7(A) = J2 XrnAiZ)a„,,z{X) . (49) 

m,Z 

This obviously depends on the composition, but is rel- 
atively insensitive to the state of ionization of the cool 
shell, as long as K-shell opacity dominates. For solar com- 
position cr(A) « 2.2 X 10-25 A^/^ cni2 HCF94|I . where the 
wavelength A is given in A. The column density of cool 
gas along the line of sight depends on the swept up mass 
behind the reverse shock according to Eqs. Q and so 
that 

7V,,„i « ^'^0 ^ - 4) M 1 

^ 1.8 X lO^sfcili.y^-if-i . (50) 

MA 

Clumping of the cool, dense shell may, however, influ- 
ence the column density (Chevalier & Blondin [l995|l . From 
their hydrodynamic calculations Chevalier & Blondin find 
a variation of A'cooi by a factor ~ 10 in different directions. 



The total emission will then be a superposition of emis- 
sion along lines of sight with different column densities. 
Since we generally do not know the degree of dumpiness, 
we, here, choose to treat the column density as a free pa- 
rameter, adjusting the absorption to flt the observations. 
For illustrative purposes we show in Fig. the effects 
of absorption using three different values for the column 
density, 10^^ — 10^^ cm^^. It is obvious that the low en- 
ergy cutoff is very sensitive to this parameter, and for a 
detailed modeling a superposition of several absorption 
components, as given by the hydrodynamic calculations, 
may be needed. 




1 10 100 

Wavelength (A) 



Fig. 14. The uppermost spectrum shows the energy emit- 
ted by the reverse shock for the same model as in Fig. 1131 
while the lower three lines represent the emerging lumi- 
nosity after absorption with column densities of iVcooi = 
10^^ cm"-^, 10^2 cm-2 and 10^^ cm-^, respectively. 



7.9. X-rays from the circumstellar shock 

While this paper mainly concerns a detailed modeling of 
the X-rays from the reverse shock, there is also a contri- 
bution from the outer, circumstellar shock. Because of the 
high velocity, and consequently high shock temperature, 
the emission from this will be dominated by continuum 
emission. This has been discussed in detail in IFLC961 A 
complication is here that the coupling between electrons 
and protons by Coulomb collisions decreases with increas- 
ing temperature, 

^'-'^ \\{)^k) Vl08cm-3J 

« 0.21 (yI^) A^Wltl days. (51) 

The electron temperature may, therefore, be consid- 
erably lower than the ion temperature. In FLC96 it was 
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found that for the parameters of SN 1993 J with a shock 
velocity of 2 x 10'' km s~^ and = 5 the ion temper- 
ature was (1 — 2) X 10^° K, while the electron temperature 
was only ~ 1 x 10^ K at 10 days. As the density is decreas- 
ing with decreasing and increasing time, the difference 
between the ion and electron temperatures can be even 
larger. A caveat here is, however, that plasma instabilities 
may lead to more efficient equipartition. The observations 
of SN 1993J, however, indicated that this was not very 
important (FLC96). For the reverse shock the density is 
so high and the temperature so low that the electrons and 
ions will generally be in equipartition. 

The spectrum and luminosity of the circumstellar 
shock depend on the mass loss rate, A^, shock velocity, 
Vs, and time since explosion, t. In contrast to the reverse 
shock, the circumstellar shock is in general adiabatic. The 
temperature behind the shock is therefore fairly uniform 
(see e.g.. Fig. 5 in FLC96). To estimate the monochro- 
matic luminosity in the X-ray band we use a similar 
method to that in IFLC96I with s = 2. Approximating 
the free-free Gaunt factor, gs, (including relativistic ef- 
fects and both ion-electron and electron-electron emission) 
with gs « 2.3T0-42A0-23 for 1 J$ A ^ 10 A, we get 



37a 22t/-1 / / 



9.6 x lO^'CAiV^ 



0.2 



^2n-0.08^-1.77g-0.144/ATg 



t 



100 days 



erg s 



'A- 



(52) 



Here Tg = Te/10^ K, A is in A, C = [H-2n(He)/n(H)]/[l + 
4n(He)/rt(H)] and / is the relative thickness of the circum- 
stellar shock region, / ~ {Res — i?s)/i?s • For ry = 7 we have 
/ « 0.3, and for ^ 12, / « 0.22 (CF94|l . The luminosity 
from Eq. H52(l should be added to the flux from the reverse 
shock, reduced by the factor below, in Eq. (|53(l . 

The ejecta are in general optically thick, and absorb 
the emission from that part of the interaction region which 
is occulted by the ejecta. For the reverse shock, where 
the cooling region is thin, the difference in radius is so 
small that we can approximate the observed luminosity 
as one half of the emitted luminosity. This has been taken 
into account when plotting the model spectra. For the 
outer shock, on the other hand, the size of the shocked 
region can be as much as 20 — 30% of the total radius, and 
this approximation is inaccurate. The relation between the 
observed and emitted luminosity can then be expressed as 



^A.obs 



0.5 + 0.5 



((/ + 1)2-1)-V2 



(53) 



For / — 0.25 this gives iobs = 0.72Lom- Thus, in 
Fig. 1151 we have reduced the outer shock emission by this 
factor, while the reverse shock emission has been reduced 
by one half, as in all other spectra. 

In addition to the free-free emission, inverse Compton 
scattering of the photospheric photons by the thermal. 



semi-relativistic electrons behind the circumstellar shock 
may contribute to the X-ray emission (Fransson I1982L 
Lundqvist & Fransson H^M FLC96). The flux of this 
depends sensitively on the temperature of the electrons, 
the mass loss rate, and the optical and UV flux of the 
supernova. In general, this component is only important 
for electron temperatures ^ 10^ K, and during the first 
months after explosion, and then only for high mass loss 
rates. This was discussed in some detail for SN 1993J 
(FLC96), where it was concluded that because of the slow 
equipartition between electrons and ions, the electron tem- 
perature was too low to scatter the photospheric photons 
up to the X-ray regime. A characteristic signature of a 
Compton component would be a powerlaw component to 
the X-rays. 

Also inverse Compton by the relativistic electrons, re- 
sponsible for the radio emission, may contribute to the 
X-rays. This component is, however, only likely to be im- 
portant for low values of A* (Chevalier et al. 12005(1 . 

As an example we show in Fig. 1151 a full spectrum 
including the thermal emission from both shocks, and the 
absorption by the shell, for A, = 1, Tg = 1, keV, V4 — 1, 
rj 13.0 and iVcooi = 10^^ cm'^ at 100 days. At this 
epoch the ion temperature behind the circumstellar shock 
is ~ 2 X 10^ K, but the electron temperature is only ~ 
(1 — 4) X 10^ K, due to the inefficient Coulomb coupling at 
the low density, ~ 2 x 10^ cm~'^, behind the circumstellar 
shock. The cutoff in the hard spectrum is therefore at 
~ 0.4 A, or - 30 keV. 

For A ^ 3 A (E ^ 3 keV) the spectrum is dominated 
by the circumstellar shock. This is also the case above the 
cutoff wavelength from the cool, dense shell at ~ 20 A. The 
figure illustrates the importance of taking both the absorp- 
tion of the cool, dense shell and the circumstellar compo- 
nent into account when analyzing the spectrum. Note also 
that the column density of the cool shell is quite uncer- 
tain fSect. [7!H|) . In addition, absorption from the interstel- 
lar medium, and possibly also the circumstellar medium, 
should be included. 



8. Conclusions 

This paper discusses for the first time a detailed, self- 
consistent model for the X-ray emission from supernovae 
interacting with a dense circumstellar medium. The model 
which we have presented here is intended to fill a gap in 
the modelling of these supernovae, and is applicable to 
the case where the reverse shock is radiative. Our model 
differs from earlier models in that the spectral code is cou- 
pled to a hydrodynamical code, and that we solve the full 
multi-level statistical equations of each ion to calculate 
the line emission. The cooling and ionization, which are 
computed for each temperature and density, are used in 
the hydrodynamical equations to determine the size of the 
emitting region. Thus the total emission from each zone is 
determined self-consistcntly, giving the correct weight to 
the emission from each zone. 
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Wavelength (A) 

Fig. 15. Total spectrum for V4 = 1, = 1, -q = 13.0, 
To = 1 keV, t = 100 days and solar composition. The solid 
line is the sum of the emission from the circumstellar and 
reverse shocks, assuming an absorbing column density of 
A^^qqj = 10^^ cm~^ in the cool shell. The dotted Hue is 
the unabsorbed contribution from the reverse shock, and 
the dashed line is the contribution from the circumstellar 
shock. 



We have shown that for a sufficiently steep ejecta den- 
sity gradient the radiative model can be applicable up to 
several years, depending on the progenitor's mass loss rate. 
We have also shown that its applicability depends strongly 
on the composition of the cooling gas, and that the high 
metallicity shocks stay radiative much longer than a shock 
with a solar composition. 

A stationary shock has been found to be a good ap- 
proximation for the cases of interest to us, while the ioniza- 
tion must be considered in a time-dependent formulation 
to account for the non-equilibrium recombination in the 
cooling gas. This effect is, however, mainly important for 
lines with A <; 100 A. We have shown that the emerging 
spectrum is highly sensitive to both the temperature of the 
reverse shock and the composition of the gas. Because of 
the strong dependence of the emerging spectrum on the 
composition of the gas, the chemical abundances in the 
ejecta can be determined from the spectrum. A require- 
ment for this is, however, a self-consistent modeling of the 
hydrodynamic and radiative processes, as demonstrated 
in this paper. 
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Appendix A: Atomic data 

The elements included in the model are H, He, C, N, O, 
Ne, Mg, Si, S, Ar, Ca, Fe and Ni. Table EU lists the 
ions for which level population calculations are performed, 
along with the number of levels included for each ion. 

The ionization rates (coUisional and autoionization) 
are computed according to Arnaud & Rothenflug H1985|l . 
except for Fe, where we use the method suggested by 
Arnaud k Raymond H1992fl . 

The recombination rates are evaluated as the sum of 
the radiative and dielectronic recombination coefficients, 
a — ar -\- ad- For the radiative recombination rates of 
C, N and O we use the fits tabulated by Nahar l(TU^ 
HM^ and Nahar & Pradhan (TWT), except for C-like N 
and O (see below). For C-like, 0-like, F-like and Ne-like 
ions the recombination rates are evaluated from the fitting 
parameters of Zatsarinny et al. (201731 I^U4all2004b|l . using 
the fitting formula suggested by Verner & Ferland H1996|l : 

ar{T) = a[(T/To)i/2(l + (7/70)1/2)1-" 

X (l + (r/ri)i/2)i+6]-i. (A.l) 

The same formula is used for H-like through Ne-like ions 
of Mg, Si, S, Ar, Ca, Fe, and Ni (Gu HUU^ and for H- 
like. He-like, Li-like and Na-like ions of the remaining el- 
ements (Verner & Ferland 11996(1 . For all remaining ions 
we use the rates compiled by D. Verner 2, which include 
data from Aldrovandi & Pequignot l|1973|l . Shull & Van 
Steenberg H1982|l . Arnaud k Raymond H1992|l . Landini & 
Monsignori Fossi (WM HW^i . Pequignot et al. lfTMT|l and 
Verner & Ferland (199gjl . 

The dielectronic recombination rates are fitted with 

MT) = ^,t^.e--'''". (A.2) 

i^l 

^ Recombination rates are available from 
,http://www.pa.uky.edu/,-^verner/fortran.html. 
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For C-like, 0-like, F-like and Ne-like ions we use the fitting 
parameters of Zatsarinny et al. H2003[l2004all2004b|l ^, and 
for the remaining ions we use the data from Mazzotta et 
al. l(TM5|l . 

The coUisional and radiative rates for the excited levels 
are taken from version 4.2 of the Chianti database (Dere et 
al. 119971 Young et al. 12003(1 ^, except for the recombination 
to individual levels. We include radiative recombination to 
He-like, Li-like and Na-like ions, using data from Verner & 
Ferland 1)1996(1 . and dielectronic recombination to Fexvi 
and Fexvil, with data from Dasgupta 1(1995(1 . 



^ Dielectronic recombination rates are available from 
|http:/ /homepages. wmich.edu/ -^gorczyca/drdata/. 

* CHIANTI is a collaborative project involving the NRL 
(USA), RAL (UK), and the Universities of Florence (Italy) 
and Cambridge (UK). 
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Table 1. Fractional abundances by number and atomic weights for the different compositions. 
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Table 2. Fitting parameters for the cooling function of Eq. 1)38(1 and the thickness parameter, Dq of the cooling region 
from Eq. 



Composition 


A (er; 


5 cm'^ 




B (er; 


5 cm'^ 




a 





Do 


(cm) 


Solar 


8.0 


X 


10" 


-23 


2.3 


X 


10" 


-24 


0.90 


0.50 


1.82 


X 


10^^ 


Helium zone 


5.5 


X 


10" 


-22 


1.5 


X 


10" 


-23 


1.10 


0.40 


1.50 


X 


lO^* 


Carbon zone 


1.9 


X 


10" 


-20 


1.3 


X 


10" 


-22 


1.35 


0.40 


1.20 


X 


lO^* 


Oxygen zone 


4.5 


X 


10" 


-20 


2.8 


X 


10" 


-22 


1.00 


0.40 


4.61 


X 


10" 



Table 3. Line strengths relative to the Oviii A18.97 line for a shock model with To ~ 1.0 keV, a single temperature 
model at the same temperature and a single temperature model at the temperature, 0.7 keV, which gives the best fit 
to the spectrum of the shock model. 



Ion 


Wavelength (A) 


F/F(Oviii) 
Shock model 


F/F(Oviii) 
Single temperature (1 keV) 


F/F(0 viii) 
Single temperature (0.7 keV) 


Cvi 


33.73 


0.15 


0.10 


0.10 


Nvii 


24.78 


0.07 


0.05 


0.05 


Ovii 


21.60 


0.23 


0.02 


0.03 


Nex 


12.13 


0.38 


0.72 


0.65 


Mgxii 


8.42 


0.14 


0.60 


0.29 


Si XIII 


6.67 


0.16 


0.65 


0.34 


Sxv 


5.00 


0.05 


0.35 


0.09 


Fexxi 


12.29 


0.46 


1.60 


1.15 



Nymark et al.: Circumstellar Interaction 
Table A.l. Number of levels and transitions included in the model for each ion. 



25 



Ion 


Lev./Trans. 


Ion 


Lev./Trans. 


Ion 


Lev./Trans. 


Ion 


Lev./Trans. 


Hi 


25/75 


Mgv 


10/25 


Sxiv 


40/117 


Fex 


172/3959 


He I 


16/187 


Mg VI 


23/142 


Sxv 


49/166 


Fcxi 


47/168 


Hen 


25/68 


Mg VII 


46/389 


Sxvi 


25/68 


Fexii 


41/696 


Cii 


18/40 


Mg VIII 


125/654 


Ar IV 


30/165 


Fexiii 


27/87 


Cm 


20/53 


Mg IX 


46/313 


Ar VII 


16/35 


Fexiv 


40/227 


Civ 


15/18 


Mgx 


40/123 


Ar VIII 


21/34 


Fexv 


53/935 


Cv 


49/165 


Mg XI 


49/166 


Ar IX 


89/1175 


Fcxvi 


21/34 


Cvi 


25/68 


Mg XII 


25/68 


Arx 


3/3 


Fexvii 


89/463 


Ni 


26/72 


Sin 


15/30 


Ar XI 


10/24 


Fe XVIII 


113/176 


Nil 


23/86 


Sim 


20/67 


Ar XII 


72/1776 


Fexix 


92/1276 


Niii 


20/56 


Si IV 


21/34 


Ar XIII 


15/44 


Fexx 


86/1309 


Niv 


20/59 


Siv 


27/115 


Ar XIV 


125/654 


Fexxi 


290/6435 


Nv 


15/18 


Si VI 


3/3 


Arxv 


46/136 


Fexxii 


204/6820 


Nvi 


49/173 


Si VII 


10/25 


Arxvi 


40/139 


Fe XXIII 


58/258 


Nvii 


25/68 


Si VIII 


72/2061 


Ar XVII 


49/166 


Fc XXIV 


40/155 


On 


15/48 


Si IX 


46/403 


Ar XVI 1 1 


25/68 


Fcxxv 


49/168 


Oiii 


46/335 


Six 


125/654 


Caix 


16/35 


Fc XXVI 


25/68 


Oiv 


125/654 


Si XI 


46/134 


Cax 


21/34 


Nixii 


31/101 


Ov 


20/57 


Si XII 


40/134 


Caxi 


89/971 


Nixiii 


48/191 


Ovi 


40/107 


Si XIII 


49/166 


Caxii 


3/3 


Nixv 


27/87 


Ovii 


49/173 


Si XIV 


25/68 


Caxiii 


20/63 


Nixvi 


40/246 


VIII 


25/68 


Sii 


43/284 


Caxiv 


15/96 


Ni XVII 


16/35 


Nell 


3/3 


Siii 


49/366 


Caxv 


46/421 


Ni XVIII 


21/34 


Neiii 


58/1212 


SiV 


52/220 


Caxvi 


125/632 


Nixix 


89/336 


Nciv 


22/91 


Sv 


16/35 


Caxvii 


46/136 


Nixx 


113/176 


Nev 


49/408 


Svi 


21/34 


Caxviii 


40/136 


Nixxi 


58/1265 


Nevi 


125/550 


Svii 


27/111 


Caxix 


49/166 


Nixxii 


15/96 


Ne VII 


46/127 


Sviii 


3/3 


Caxx 


25/68 


Ni XXIII 


20/79 


Ne VIII 


40/103 


Six 


86/2640 


Fell 


142/1268 


Ni XXIV 


125/654 


Ne IX 


49/173 


Sx 


22/215 


Fc VI 


80/2083 


Nixxv 


46/136 


Nex 


25/68 


Sxi 


46/416 


Fe VII 


9/23 


Ni XXVI 


40/161 


Mgii 


21/34 


Sxii 


125/654 


Fe VIII 


83/267 


Ni XXVII 


49/168 


Mg IV 


3/3 


Sxiii 


46/136 


Feix 


13/45 


Ni XXVIII 


25/68 



